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We report the first observation of stable single photon sources in silicon carbide 
(SiC). These sources are extremely bright and operate at room temperature demon- 
strating that SiC is a viable material in which to realize various quantum information, 
computation and photonic applications. The maximum single photon count rate de- 
tected is 700k counts/s with an inferred quantum efficiency around 70%. The single 
photon sources are due to intrinsic deep level defects constituted of carbon antisite- 
vacancy pairs. These are shown to be formed controllably by electron irradiation. The 
variability of the temporal kinetics of these single defects is investigated in detail. 




Silicon carbide has fast become an important 
wide-bandgap semiconductor for high-power, high- 
temperature applications and is regarded as a key ma- 
terial for next generation photonic^ and electronic 
devices^'^. As such, SiC has well-established growth and 
device engineering protocols^. This has allowed the de- 
velopment of a mature platform for nano-fabrication and 
high-Q photonic crystal cavities, which have already been 
reported for a spectral region from 550 to 1450 nm^. SiC 
is considered one of the best biocompatible materials, es- 
pecially to blood^, and the synthesis of SiC nano-particles 
have been employed as fluorescent biological labels in cell 
imaging^. 

Such versatility has led to recent interest in the pos- 
sibility of using SiC as a platform for various photonic, 
spintronic and quantum information applications^. In 
particular, it is a promising material in which to host 
qubits and thus holds great promise for the future imple- 
mentation of scalable quantum processors^' In partic- 
ular, several authors^'^^'^^ indicate that SiC could harbor 
equivalent defects to the nitrogen vacancy center in dia- 
mond, a well known solid-state qubit^^. While candidate 
ensemble optically- active defects, having the net angular 



momentum (or spin) required to act as qubits operating 
at room temperature, have been identified^ single 
photon sources and the ability to engineer them control- 
lably have not yet been demonstrated in SiC. 

In this work, we present, to the best of our knowledge, 
the first observation and engineering of single photon 
sources from point defects in SiC. At room temperature 
the defect responsible emits between 670-710 nm and is 
assigned to the carbon antisite- vacancy pair^^'^^. A great 
number of these single defects can be created, indicat- 
ing that the electron irradiation results in a very efficient 
method to engineer single defects in SiC. The study of the 
kinetics and photoluminescence (PL) of these single de- 
fects using confocal microscopy and single photon corre- 
lation, reveals a variability of the photophysics of the de- 
fects, that is undetectable with ensemble measurements. 
As the identification of deep defects in wide bandgap 
semiconductors is typically carried out by photolumi- 
nescence and electron paramagnetic resonance ensemble 
measurements^^, the direct correlation to single defects 
is often difficult. Here, we perform a study of the defect 
formation from the ensemble level down to the single de- 
fect level, providing a reliable route for the identification 
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of the defect responsible for the observed single photon 
source. We also correlate the carbon antisite- vacancy 
pair to the formation of another interesting defect for 
quantum information, the silicon vacancy in SiC^^. 

The discovery of single photon emission in this ma- 
terial and its efficient creation, is the starting point to 
pioneer a new branch of alternative defects for next gen- 
eration solid state candidates in quantum information, 
computation and photonics. 



RESULTS 

Single defect measurements 

Figure la,b shows confocal PL maps of a high purity, 
semi-insulating SiC wafer irradiated with 2 MeV elec- 
trons and subsequently annealed (see Methods: Samples 
preparation). Spatially isolated, extremely bright defect 
centers are observed, the number and brightness of which 
are dependent on the excitation wavelength as shown in 
the histogram of Fig. Ic. It is clear that the 532 nm ex- 
citation laser allowed a greater variety of defects to be 
addressed. As observed in the confocal maps in Fig. 1, 
higher count rates are detected when using the 532 nm 
excitation laser, mainly due to the use of an ultra-steep 
long pass filter with 90% transmittance so that the sin- 
gle photons were collected on their entire spectral band- 
width. Centers could also be observed in untreated sam- 
ple but at a concentration of almost an order of magni- 
tude lower than the processed sample (Fig. Ic). The con- 
centration of these centers was found to increase mono- 
tonically with electron irradiation fluence and increas- 
ing anneal temperature (see Supplementary). Using a 
Hanbury-Brown and Twiss interferometer to measure the 
photon correlation we found that all the diffraction lim- 
ited spots with diameters of around 300-400 nm demon- 
strated single photon emission with wavelengths between 
670-710 nm. 

Typical PL spectra of these bright centers are shown 
in Fig. Id excited using both 532 and 660 nm lasers. The 
width of the PL emission tended to be narrower when 
excited with a 660 nm laser. However, the PL line-shape 
possesses some variability both at 532 nm and 660 nm ex- 
citation as shown in the Supplementary Information. Re- 
markably, photon correlation spectroscopy revealed that 
centers emitting single photons were the dominant de- 
fect, i.e. the second order photon correlation function at 
zero delay is mostly ^^^^(0) < 0.5 (Fig. le). 

The photo-physics and the population dynamics from 
the ground state to the excited state of these single pho- 
ton sources was investigated. The PL temporal dynamics 
can be studied by analyzing the g^'^\r) at different ex- 
citation optical powers. Fig. 2a shows the anti-bunching 
behavior of the center with PL shown in Fig ld(i) for a 
range of 532 nm excitation powers. Due to the presence 
of a photon-bunching effect (higher correlation at longer 
delay times), each curve was fit with a three- level system 



model, 

g^^\r) = 1 - (1 + a)exp[-Air] +aexp[-A2r] (1) 

where the decay rates, Ai and A2 are free parameters 
and related to the transition rates, r^j, from level {i) to 
level (j) by the equations Ai = ri2 + r2i, A2 = r^i + 
^23^12 /Ai^^. The values of a, is a measure of the photon- 
bunching, that prevents radiative decay to the ground 
state (1) from the excited state (2) by a metastable state 
(3) and can be expressed as a = ri2r23/(Air3i). Photon- 
bunching at longer delay times in single molecules^^ and 
color centers^^ in diamond has been associated with the 
presence of such a long-lived metastable state. Non- 
radiative decay to the ground-state takes place from this 
metastable state, inducing an overall reduction of the sin- 
gle photon emission arising from low rate electron trap- 
ping in the metastable state. 

For this particular center the saturation count rate and 
optical saturation power are = 1.3x 10^ counts/s and 
Psat = 0.5 mW, respectively. This is obtained from fit- 
ting the data in Fig. 2b with <l> = ^ooPopt/{Psat + Popt)- 
The dependence of the transition rate parameters of this 
model, a and Ai^2, on the excitation power are shown 
in Fig. 2(b,c). A linear optical power dependence of 
the coefficient Ai = r2i(l + o^Popt) is assumed, where 
ri2 = ar2iPopt = ^^opt/ ihv)^ is related to the absorp- 
tion cross- section of the defect, a. Due to the very weak 
dependence of the metastable state population on the ex- 
citation power, we have assumed that the shelving (r23) 
and de-shelving (rsi) decay rates of the metastable state 
are constant. Other models used to describe the popula- 
tion dynamics of the excited state and of the metastable 
state versus the excitation optical power^^'^^, however, 
may also describe our data. For example, in some cases 
a linear dependence of rsi = r2,i(l + PPopt) may be em- 
ployed where /3 is a fitting parameter that describes the 
variation of the decay rate from the metastable to ground 
states with changes in optical excitation power. This pa- 
rameter is important since it signifies that a power de- 
pendent effect in de-shelving the metastable level is oper- 
ative. This dependence has also been observed for single 
molecules^^. Fits with and without the f3 term are shown 
in Fig. 2(b,c). 

On fitting the three-level system model, we obtain the 
decay rate terms r2i, r^i^ r23, a and /3. Extrapolation to 
zero excitation power gives the excited state lifetime, r2i, 
and the shelving decay rates, rai. While extrapolation to 
infinite excitation power of A2 gives r23. For the center 
shown in Fig. 2 we obtain l/r2i = (1.2 ± 0.1) ns, a = 
(0.3 ± 0.1), /3 = (2.1 ± 0.5) mW-\ = (45 ± 3) ns, 

1 A23 = (7 ± 2) ns, and for /3 = 0, l/rg^ = (23 ± 2) ns, 
1/^23 — ^ 2) ns. The absorption cross section is a ~ 
5 X 10~^^ cm^. Using these values the following model 
can be used to determine the quantum efficiency, r]qe, of 
this single photon source, 

^ = Vqer]collVdetr2l{l + ^21/^2 + r23/r3i)"^ (2) 
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FIG. 1. Confocal microscopy images and spectra of single photon sources in SiC. Confocal map of the low electron 
fluence sample excited with 660 nm (a) and 532 nm (b) lasers. The diffraction limited spots all correspond to single photon 
emitters. Map (a) has been scanned using a 10 nm FWHM 694 nm filter. Note the difference in color scale, c, Histogram 
of the number of centers observed as a function of the spatial FWHM measured and counted from (a,b). The histogram 
from a confocal map of an unprocessed sample is also shown indicating a clear increase in the defect density after irradiation 
and annealing at 300° C. d, Normalized room temperature PL of a single photon source in the irradiated sample excited with 
532 nm (i) and 660 nm (ii,iii) lasers. The laser wavelengths are indicated by the dashed lines, d(iii), PL of a center which was 
not anti-bunching. The transverse and longitudinal optic Raman modes (TO, LO) are also indicated, e, The corresponding 
anti-bunching curves for d(i) and (ii). 



By fitting the count rate we obtain rjqeTjcoii^det = 0.008 ± 
0.001. If we assume that the detection efficiency is only 
limited by the detector's quantum efficiency (/?det=0.55) 
with the filters having a transmission of 90%, and a col- 
lection efficiency of our setup determined to be rjcoii ^ 
0.02%^'^, r]qe ^ 70% is estimated. 

To assess the properties of the single photon emis- 
sion accurately, the variability of their kinetics must 
be established. Therefore, many photo-stable emitters 
were measured with the two different excitation wave- 
lengths and analyzed. We observed a moderate disper- 
sion of the excited state lifetime from 1.2-2.1 ns, yield- 
ing average values of 1.5 ± 0.3 ns (20% variability) and 
a = (0.4 ± 0.3) mW~^. Concerning the metastable 



state, due to the very small photon-bunching effect ob- 
served with increased excitation power, the measurement 
of the decay rates from and to the metastable state are 
affected by a large variability, giving average values of 
1/rgi = (40 ±20) ns and l/r23 = (13±5) ns. In addition 
not all the emitters were best fit with the (3 parameter. 

The saturation properties of the defects also display 
a remarkable variability, with optical saturation powers 
observed between 0.2 to 3.3 mW and saturation count 
rates found from 2.5 x 10^ count s/s to the maximum 
observed of 2 x 10^ counts/s. 

The observed SiC defects have a much shorter 
metastable state lifetime compared to typical diamond 
color centers (of the order of 100 ns for NV or jus for 
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FIG. 2. Temporal dynamics of single photon sources in SiC. a, Second-order correlation function at different 532 nm 
excitation powers from the center with PL shown in Fig. ld(i). A weak dependence of the metastable state shelving with the 
optical excitation power is observed. The curves are shifted vertically for clarity, b, Saturation curve of the defect fitted with 
Eq. 2 where /3 = 2.2 mW~^ (solid line) and /3 = (dashed line), to infer the quantum yield (left axis), a versus excitation 
power, fitted with the /S term (solid) line and with /3 = (dotted line) (right axis), c, Decay rates, Ai and A2 extracted from 
the fits in a versus the excitation power fit with the model with the /3 term (solid line) and without it (dashed line). 



other defects such as Cr-related or SiV centers) In 
addition, the here observed large dispersion of its lifetime 
could indicate that the photon-bunching effect in the cor- 
relation function may be also attributed to a longer-lived 
non-radiative state with different charge or to thermal 
excitation from the excited states to higher energy non- 
radiative states. These states can either belong to the 
center or they can be formed by other impurities in the 
local nanoenvironment, as it occurs in some diamond de- 
fects. Low temperature and longer delay time photon 
correlation measurements are required to fully identify 
the photon-bunching associated physics. 



Single photon source stability 

We studied the PL emission stability of these defects 
and observed that 45% of the emitters measured in the 
irradiated and 300° C annealed sample are stable over 
a long period of time. An example of a center with 
a stable PL time trace is shown in Fig. 3a along with 
its associated Gaussian count rate histogram Fig. 3(b). 
Interestingly, the remaining centers show PL blinking 
with some eventually photo-bleaching, even if not perma- 
nently. Such centers in samples annealed at higher tem- 
peratures (500° C) were almost absent suggesting that the 
blinking was related to the density of non-radiative im- 
purities which could be annihilated with annealing. An 
emitter with on-off fluorescence statistics, between two 
PL states is shown in Fig. 3(c). By measuring all on/off- 
time durations in a single trace, the corresponding prob- 
ability of occurrence P{ron,off) can be calculated^^. The 
distribution of the "on" and "off" times are plotted on a 
semi-logarithmic scale in Fig. 3(e) together with the fits. 
Both Pijon) and P{roff) show an exponential trend fol- 



lowing the equation: 

P{ron,off) oc exp[-r/(ron,o//)] (3) 

with the time constant. Ton usually much greater than 
To//, especially in unprocessed SiC or samples annealed 
at higher temperatures. This model suggests that typi- 
cally a single mechanism is operative, which could involve 
a photo-induced charge conversion of the addressed cen- 
ter due to defects in close proximity acting as charge traps 
or donors. P{Ton,off) is also found to have a strong laser 
power dependence supporting a photo-induced mecha- 
nism (see Supplementary). Two state intermittence is 
the most common blinking observed. However, a small 
number of centers displayed more complex PL blinking. 
This was especially true when using 532 nm laser ex- 
citation where three states or emission with no stable 
level could be observed from a single anti-bunching cen- 
ter. Centers excited with the 660 nm laser showed greater 
photo-stability 

Ensemble measurements 

Ensemble room temperature and low temperature 
measurements were performed on a sample irradiated to 
10^^ e/cm^ at 532 nm using a Renishaw Micro- Raman 
spectrometer. At 80 K the so called AB lines can clearly 
be observed (Fig. 4b, c) and result in a broad PL band 
at room temperature. These lines are associated with 
the neutral (CsiVc)^ defect The atomic model of this 
defect with its four possible cubic or hexagonal orienta- 
tions is shown in Fig. 4a. The single photon emissions 
observed in this work are attributed to this particular de- 
fect. The 532 nm laser excites all zero phonon lines. The 
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FIG. 3. Photo-stability of single photon sources in SiC. a, 20 s fragment of a PL time-trace sampled into 10 ms bins 
from a stable single photon source with its corresponding count rate histogram (b). c,d, A similar time trace is shown for a 
blinking center. The on and off times of (b) on a semi- logarithmic scale, f, Histogram of the measured Ton and To// of the 55 
emitters analyzed. 45% of the emitters in this sample were stable and did not show blinking behavior. 



660 nm laser excites only the B spectra and the phonon 
side band is detected. 

Similarly, the VI and V2 lines^^ are also observed, cor- 
responding to defects on inequivalent sites of the Vsi de- 
fect (Fig. 4b, d). These two defect types, CsiVc and Vsi, 
are the dominant defects observed in this study. Confo- 
cal mapping revealed that the VI and V2 lines are domi- 
nant and have a higher intensity than the AB lines in the 
high electron fluence samples. In samples irradiated with 
lower fluences this defect was rarely observed due to the 
greater threshold energy required for the Si sub-lattice 
of SiC. We did not observe anti-bunching in the Vsi cen- 
ter possibly due to the dominance of the CsiVc defect. 
Excitation with wavelengths closer to the V1/V2 system 
may be required to observe single photon emission even 
in un-processed samples. 



DISCUSSION 

These experiments have demonstrated the formation 
of stable single photon sources in SiC operating at room 
temperature. The emission around 700 nm can be mod- 
eled with a three-level system containing a metastable 
state with a short lifetime compared to single photon 
sources in diamond. The brightness of this single pho- 
ton source is comparable or even higher than the bright- 
est color centers observed in bulk diamond^ ^. Whilst 
as grown defects yielding single photon emission are ob- 
served, scaling of the number of defects with irradiation 
fluence, indicate that the defects can be successfully en- 
gineered. Unusual blinking behaviour is observed which 
is possibly due to non-radiative recombination at other 
defects sites in close-proximity to the addressed center. 
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FIG. 4. Ensemble measurements of defects in electron irradiated SiC. a, Structure of the four possible CsiVc type 
defects in SiC shown in the (1120) crystallographic plane. Small and large circles represent C and Si atoms while the small open 
circles are vacancies, h and k denote the hexagonal and cubic sites in the lattice, respectively, b, PL of a sample irradiated 
with 2 MeV electrons to a fluence of Ix 10^^ cm~^ at temperature of 80 K (top spectrum) to room temperature (bottom 
spectrum). Expanded view of the 80 K PL spectra showing the CsiVc (c) and Vsi (d) defect signatures. The asterisk in (c) 
indicates a signal of unknown origin. 



Blinking could be reduced when a high annealing tem- 
perature was employed. The defect responsible for the 
single defect emission is identified as the carbon anti- 
site vacancy pairs in its neutral charge state, (CsiVc)^. 
The deployment of these single centers and their further 
characterization is extremely promising in quantum in- 
formation applications and may trigger further studies of 
point defects in SiC at the single defect level. 



METHODS 

Sample preparation 

The material under study is high purity semi- 
insulating (HPSI) SiC purchased from CREE. According 
to the manufacturer, the as-received samples contained 
B (IQi"^ cm-3), C interstitials (10^^^ cm-^) and N (<10^'^ 
cm~^). 

Samples were irradiated with 2 MeV electrons to flu- 
ences of 1x10^^, 1x10^^ and Ix^'^ e/cm^ at a flux of 
2.82x10^^ cm~^s~^. The samples were held at tempera- 



tures below 80° C during irradiation in a nitrogen ambient 
to inhibit Frenkel pair annihilation and reduce the pos- 
sibility of contamination. After irradiation, a low tem- 
perature (300° C) anneal in an Ar ambient for 30 min- 
utes was performed. This temperature was chosen so 
as to enhance the concentration of various defects of in- 
terest based on ensemble electron spin resonance (ESR) 
measurements^^. An unprocessed sample from the same 
wafer was also characterized in order to establish the ef- 
fectiveness of the irradiation procedure employed. 

Similar samples were also annealed for 30 minutes at 
300, 400, 500° C in an Ar ambient after 2 MeV electron 
irradiation to a fluence of 10^^ e/cm^. 



Photoluminescence and photon correlation 
measurements. 

A Hanbury-Brown-Twiss interferometer was used to 
identify optically active single defects and measure the 
time correlation of the center emission. In the confocal 
microscope the samples were excited using 0.2-1.5 mW 
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continuous wave laser operating at 532 nm (Coherent, 
model: Compass 315M-100) and continuous diode laser 
operating at 665 nm. The laser was focused onto the 
sample using a lOOX infinity-corrected oil immersion ob- 
jective lens with a NA of 1.4 and the luminescence was 
collected confocally through a 50 /im pinhole. A spec- 
trometer (Princeton Instruments Acton 25001, Acton) 
with a cooled CCD (Camera Pixis 100 model: 7515- 
0001 Princeton Instruments) was used to characterize 
the luminescence and a HBT interferometer with single- 
photon-sensitive avalanche photodiodes (Perkin Elmer 
SPCM-AQR 14) was used to measure the photon statis- 
tics. Photon counting and correlation was carried out 
using a time-correlated single-photon-counting (TCSPC) 
module (PicoHarp 300, PicoQuant GmbH). 

The SiC samples were mounted on a piezo XYZ stage 



allowing 100x100 fim^ scans. The unwanted residual 
laser line was eliminated by a 660 nm long pass filter. 
A filter with 694±10 nm was placed in front of HBT 
interferometer to observe single photons during 660 nm 
laser excitation. The spectroscopy was done using only a 
long pass filter. All the single photon/single defect mea- 
surements were performed at room temperature. 

For ensemble PL measurements we used a Renishaw 
MicroRaman Spectrometer equipped with a Linkam 
stage for temperature control. The excitation source was 
a 532 nm laser with 60 mW of power at the back of the 
objective. The objective was a 50x long working distance 
with 0.5 NA. A 532 nm edge filter was used. The spec- 
trometer has a working distance of 0.25m and uses an 
8/im slit and a 1200 g/mm grating. The detector is a fan 
cooled Si CCD. We used 10 s accumulations per point. 
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PHOTOLUMINESCENCE VARIABILITY 



We have used single defect spectroscopy to study the 
photo-physical parameters of single defects in SiC not 
accessible with conventional ensemble techniques such as 
standard electron paramagnetic resonance. A detailed 
analysis of the variability of these parameters is made 
possible. Firstly, the photoluminescence (PL) spectra of 
single photon sources in SiC is shown in Fig. SI. The PL 
of samples irradiated with 2 MeV electrons to a fluence 
of 10^^ e/cm^ and annealed at 300° C are shown excited 
with either a 660 nm laser or a 532 nm laser are shown 
in Fig. SI (a) and (b), respectively. Broad PL is observed 
between 600-800 nm. Fig. SI (c) shows typical spec- 
tra from a sample irradiated to a fluence of 10^^ e/cm^ 
followed by a similar anneal. The emission clearly corre- 
sponds with the Vsi-related emission shown in Fig. 4 of 
the main article. The emission related to the (Cg^VQ)^ 
shows much greater temperature quenching but the PL 
band observed around 600-800 nm is associated with this 
defect and is responsible for the the single photon emis- 
sion. There are in fact six different zero phonon lines 
associated with this defect in SiC which contribute to 
the variability of the band peak position. 



BLINKING STATISTICS 

We studied the statistical behavior of the "on" and 
"off" time distributions of single emitters and identified 
a number of blinking categories. One common case is an 
unclear blinking behavior as shown in Fig. S2(a). Emit- 
ters having two on-off states are shown in Fig. S2(b) and 
(c), where in case (b) the "on" time dominates the "off" 
time. The opposite is true in (c) where the "off" state 
has a non-zero intensity. Finally, Fig. S2(d) shows blink- 
ing between 3 separate states. This last example is less 
common but is more likely to occur when 532 nm laser 
is used for excitation. Fig. S2(e) shows the dependence 
of the "on" and "off" times on the excitation power for a 
two state single blinking centre, both the "on" and "off" 
times reduce with increasing excitation power. 
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FIG. SI. Variability observed in PL emission at room 
temperature using confocal microscopy. PL emission 
associated with the single photon center excited with 660 nm 
(a) or 532 nm (b) in samples irradiated with 10^^ e/cm^. (c) 
PL emission measured in the sample irradiated to a fluence 
of 10^^ e/cm^, excited at 660 nm. These emissions do not 
correspond to single defects. 
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FIG. S2. Example of blinking statistics of single defects, a-d show examples of blinking of single photon emitters in a 
SiC sample irradiated with 2 MeV electrons to a fluence of 10^"^ e/cm^ and annealed at 300° C. The time constants for each 
state are also given, e, The on/off time constant dependence excitation power. 



DEFECT CONCENTRATION VERSUS 
IRRADIATION FLUENCE AND ANNEALING 
TEMPERATURE 

Using confocal imaging, we studied the concentration 
of the defects in samples prepared with different electron 
irradiation fluences 10^^, 10^^ and 10^^ e/cm^. For this 
purpose we used an excitation at 660 nm. We compared 
the results obtained with the PL observation to make 
sure we were observing the same defects. In addition, 
we also imaged a sample of the same SiC purity and 
transparency to determine preexisting as grown single 
defects. As expected the material contains the defects as 
grown with single photon emission signature, however we 
show that the number of defects scale with the electron 
irradiation fluence, indicating that the defect responsible 
for the single photon emission can be induced from the 
electron irradiation. 



We measured the concentration of the defects in sam- 
ples prepared with the same electron irradiation fluence 
of 10^^ e/cm^ , applying different annealing temperatures 
(as described in Methods). For this purpose we again ex- 
cited at 660 nm. We checked that the bright spots in the 
confocal were corresponding to single emitters with the 
same obtained PL emission. We observed a remarkable 
increase in the number of defects for higher temperature 
annealing. A summary of the total number of defects in 
a 100 X 100 /im^ scan versus temperature and fluences is 
shown in the Fig.S3(h). 

The increase in the number of defects with the in- 
creased annealing temperature is indicating that the de- 
fects studied here are more stable after high annealing 
temperature. In the 500 °C annealing sample some de- 
fects were also appearing to be more than 1 (^^^^(0) ~ 
0.5), and very few defects were blinking, even if photo- 
bleaching was still observed. 
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FIG. S3. Confocal maps with defects concentration versus irradiation fluence and annealing temperature, a- 

g Confocal maps of SiC with the processing conditions indicated. Note that only g is on a different color scale since it is 
relatively bright and shows no isolated defect PL. h The number of diffraction limited optical centers as a function of the 
anneal temperature (top axis) and fluence (bottom axis). 



